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Asymptotic expansion of the expected spectral 
measure of Wigner matrices 
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Abstract 

We compute an asymptotic expansion with precision 1/n of the moments of the 
expected empirical spectral measure of Wigner matrices of size n with independent 
centered entries. We interpret this expansion as the moments of the addition of the semi¬ 
circle law and 1/n times an explicit signed measured with null total mass. This signed 
measure depends only on the second and fourth moments of the entries. 
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1 Introduction and main result 


Let n > 1 be an integer. In this note, we will be interested in the so called Wigner ensemble of 
n x n matrices of the form 


x(") = (xg } ) 


1 <i,j<n 


X d' } ■sjncr 2 ™ 1 ’’ 


(1) 


where (W, y) ;>/>1 are independent random variables with mean 0 and moments of all order. 
We will consider two cases depending on whether the random variables are real or complex: 


Real case: we set r = 1 and the random variables are such that Wy, = Wfj and 
E[W, v ]=0, E[W?-]= X 2 

1 1 S if 1 = ] 1 


( 2 ) 


Complex case: we set r = 0 and the random variables are such that Wy j = W; r y and 


E[W irj \ = E[Wlj] = 0, E[|W f/ y| 2 ] = 


a 2 if i j 


if i = j 


!, E[ W irj ] = cc if i (3) 


Among the real and complex cases, one can find the Gaussian Orthogonal Ensemble 
(GOE, corresponding to r = 1, s 2 = 2, <j 2 = 1 and cc = 3) and the Gaussian Unitary Ensemble 
(GUE, corresponding to r = 0, s 2 = 1, a 2 = 1 and a = 2). 


Our main purpose is to study the spectral measure 
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}hi — 


n 


E h- 

AeSp (X<">) 


(4) 
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of the matrices X^ n \ More precisely, we compute an asymptotic expansion in n of the expected 
moments 


m k (n ) 



(5) 


of the spectral measure y n : 

Theorem. For every k > 0, as n —>• oo, one has 


m k (n) = I x k dsc(x) + ^ J x k dv(x) + o 


( 6 ) 


where sc denotes the semi-circle measure having density 

— V4-x 2 l {xe{ _ 2r2]} , 

and v is the signed measure with 0 total mass given by: 



Remark. For the GUE, the measure v is null, and for the GOE it is simply 


1 

2 




+ ^{x=-2}) 


1 {xe[- 2 , 2 ]}dx\ 

TZ\J 4 — X 2 J 


1 {xe[-2,2\}d x 
TCy /4 — X 2 


Before going through the proof of the theorem, let us make some comments. Asymptotic 
expansions of the expected spectral measure are of crucial interest for central limit theorems. 
Indeed, they are usually stated as follows for Wigner matrices: if / is smooth enough, the 
random quantity 


n 


j f(x)dy n (x) — E J f(x)dy n (x] 


converges to a Gaussian random variable as n —>• oo as proved by Lytova and Pastur in [j6]| . 
Elence it is important to know the centering term with a precision of order o(l/n). 

For the special case of the GUE, the measure v is null. This agrees with the full genus 
expansion of m k (n ) computed by Brezin Itzykson Parisi and Zuber |T] and Harer and Zagier 
|j2|. For the GOE, Ledoux 0 computed an analog of this full expansion which is also in 
agreement with our result. The full expansion of m k (n) is not known for other ensembles. 
However Johansson J3| computed the term of order 1/n for invariant matrix ensembles with 
a potential. This framework coincides with the case of matrices with independent entries 
only for the GUE and GOE. 

The first order in the asymptotic expansion of the expected Stieltjes transform of y n is 
computed by Khorunzhy Khoruzhenko and Pastur 0] for general matrices with independent 
entries. Their result is valid only at some finite distance from the real line. This analytic 
approach does not yield an asymptotic expansion for m k (n ) and is radically different from 
the combinatorical method used in our computations of these moments. 


Next section is devoted to the proof which enumerates some classes of closed paths which 
are fancy variants of the contour functions of trees. They are reminiscent of the combinatorics 
of trees and maps involved in the Harer Zagier expansion of m k (n ) for the GUE. 
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2 Proof of the theorem 

We start by writing the moments 


mjt(n) = E 


x k d}i n (x) 


l+7rrfc • 


n 


E E K, i 2 w W3 ... w Wl ] 


2 <- T ' ii.i t -l 


As it is done classically, we regroup words 44 • • • hh i n equivalence classes (we say that two 
words are in the same class if they coincide up to a permutation of the letters of the alphabet 
{1,..., ft}). The main reason to do so is that two equivalent words Z]Z 2 ... 44 and z^ ... i' k i[ 
have the same contribution in m k (n): 


E [w„, i 2 w i2 ,, 3 ... wu] = E W V W V3 ... w M 


If c is an equivalence class of words, we denote by E [W c ] the common contribution of its 
elements (i.e. the quantity E [W^W^ ... W^-J for one of its elements z'iz '2 ... 44). We will 
often see a word 44 ... 4 z \ as a closed spanning walk with k steps on a connected graph 
with vertex set {4 ,..., 4 } and edge set {{z'i, z‘ 2 }, { 4 , z’ 3 },..., { 4 , 4 }} (these graphs have no 
multiple edges but can have self-loops). Notice that the graphs and walks associated to two 
equivalent words are isomorphic. 


Fix k, v and e some integers, we denote by A (. k , v, e ) the set of all equivalence classes of 
words z'iz '2 ... 44 on the alphabet { 1 ,..., n} such that 


H4,---,4} = v and H{ z ’i/ z ’ 2 }/{« 2 , 4 }/---/{ 4 / 4 }} =e. 


If (, k , v, e) 7 kQ) satisfies E [W c ] 7 ^ 0 for some c G ^{k, v, e ), one always has 

v<e + l and e< [k/ 2 J. 


The first inequality is due to the fact that any connected graph with v vertices has at least 
v — \ edges. The second inequality comes from the fact that the random variables W hj are 
centered and therefore every edge of the graph defined by c has to appear at least twice in a 
word belonging to c. 

Now notice that every c G ^(k, v, e) is the equivalence class of exactly n(n — 1) ■■■ (n — 
v + 1) different words. This gives the expression 

1 IT/2J e +l 

m k( n )= k E Ew(w-i)•••(»-p + i) E E t w c]. (7) 

n 2(7 e—i v —\ cetf(k,v,e) 

We are now ready to compute the asymptotic expansion of w^(n) announced in the theorem. 
First, we prove that the momemts to the moments of the semi-circle law, next we prove that 
the magnitude of the next term in the expansion is 1 / n and finally we compute the term of 
order 1 In in the expansion. We identify these terms as the moments of the measure v given 
in the theorem by computing their Stieltjes transform. 
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2.1 Limits of the moments 

Non vanishing terms in 0 are such that v = /c/2, and therefore e = k/2 and k must be even. 
In that case every c G ( tf(k,k/2 + l,/c/2) is such that E[W C ] = ( cr 2 ) k/2 and 

W#{k,k/2 + l,k/2) = Cat(fc/2) 

yielding the moments of the semi-circle law for lim n —»• oo m^{n). This is basically Wigner's 
original proof of convergence to the semi-circle law. 

2.2 The leading vanishing term is of order 1/n 

Fix a > 0. In order for (0 to have a term of order 1/ n a , there must be 1 < v < |_ / J + 1 such 
that one has: 

n v _ 1 
n f+i n a 

and therefore 

a = l + v> k W 2\>{° 

2 ~ 2 L J - \l/2 

Now suppose k = 21 + 1. There is a term of order 1 /\Jn in 0 if and only if v = l + l and e = l. 
This means that the graph spanned by the word zp 2 • • • hh is a t ree with l edges. In addition, 
the word is a closed walk of length 21 + 1 on this tree, which is impossible. Therefore, the 
leading vanishing term in 0 is at least of order 1/n. 

Notice that this also means that if k is odd, the leading term in 0 is of order smaller than 
or equal to n“ 3/2 . 

2.3 The term of order 1/n 

We suppose now that k = 2Z is even and compute the contributions of order 1 /« in 

m 2 i(n) = : Y^ n (n- 1) •••(« — P + 1) E t W c] ( 8 ) 

n a v,e cetf{2l,v,e) 

where the sum on v, e is taken over v and e satisfying 

1 <v < e + 1 and 1 < e < k. 

The values of v and e yielding terms of order 1 /n in 0 are: 

1. v — e + 1 — l + l; 

2. v = e + 1 — l; 

3. v = e = l. 


if k even; 
if k odd. 


Case l:v — e + 1 — l + l. 


In this case, the contribution of order 1/n is the product of the coefficient of n v ~ 1 in the 
product n(n — 1) • • • (n — v + 1) and Cat(l ) as computed above for the order 0. The total 
contribution of these terms is therefore 



1Z(Z + 1) 

n 2 


Cat (Z). 
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Figure 1: Labeled tree with contour S. 


Case 2 : v = e+1 = l. 

In this case, the graph defined by a class c G ^ (21, v, e) is a tree with l — 1 edges. In opposition 
to the previous case and the computation of order 0 , where each edge of the tree was visited 
exactly twice, here, exactly one edge is visited four times and all the other edges are visited 
twice. Note also that because the graph is a tree, if an edge ij is visited in that order, the next 
visit of this edge must be done in the reverse order ji. 

Consequently, we have to enumerate sequences of the form 

S = i\ Si i jS 2 ] iS 3 i jS^ z'Ssh 

such that the sequences i \ Sp'S^q, /S2/, 2S3Z and /S4/ are the contour functions of disjoint trees 
with respectively p\, pi, p 3 and p 4 edges satisfying p\ + pi + po> + P 4 = k ~ 2 . 

The sequence z’i S1/S5Z1 corresponds to a rooted tree with p\ edges and a marked corner 
(adjacent to the vertex i) where the trees corresponding to the three other sequences are 
inserted (see Figure [T] for an illustration). Note that when p\ = 0 , the sequence S boils down 
to f yS2y iS^ijS^j i. 

There are (2 pi + l)Cat(pi) such trees with a marked corner. In that case every c G 
^(21,1,1 — 1) is such that E[W C ] = (cr 2 y~ 2 a. and 

§<#(21,1,1 - 1) = (2pi + l)Cat(pi) • • • Cat(p 4 ). 

PlH-hp4=/—2 

The contribution to the term of order 1/n in the asymptotic expansion of the moment mzi °f 
this case is then: 

m ii = E ( 2 Pi + !) Cat (Pi)' • ’ Cat(p 4 )- 

Pi +-+P4=J-2 


Case 3: v = e = l. 

In contrast with previous cases, there are two kinds of equivalent classes in <#(2l, 1,1) leading 
to two different contributions E[W C ]. We detail each kind in the following. 

First kind: graphs with a self-loop. 

Suppose that for some 1 < j < 21 we have ij = z’y_|_i. We have to enumerate sequences of the 
form 

z" 4 S 4 z" 1 S 211 S 3 I 1 
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and such that the sequences q SpSsq and i S 2 z are the contour functions of trees with respec¬ 
tively pi and p 2 edges satisfying p\ + p 2 = l — 1. 

For this kind of equivalence class, E[W C ] = (o' 2 ) 1 's 2 and there are 

I] (2pi + l)Cat(pi)Cat(p 2 ) 

Pl + P2=l-1 

such equivalent classes. Hence, the contribution to the term of order 1 / n in the asymptotic 
expansion of the moment 11121 of this case is: 

m ii = ““2 E (2pi + l)Cat(pi)Cat(p 2 ). 

Pi+p 2 =z-i 


Second kind: graphs with no self-loop. 

In that case, the underlying graph spanned by the walk qz 2 ■ ■ ■ hih is a tree with l — 1 edges 
and an additional edge between two of its vertices forming a cycle. The graph thus consits 
on a cycle of length p > 3 and trees attached to vertices of the cycle. We denote by C \,..., c p 
the vertices of the cycle. 

First, notice that since e = l, each edge is visited exactly twice by the walk. We distinguish 
two cases: 

1. edges of the cycle are visited one way (twice from z to /); 

2. edges of the cycle are visited both ways (once from i to j and once from j to i as in the 
previous cases). 

The first case contributes only for real matrices and corresponds to sequences of the form 

h Li C\C 2 L 2 C 2 C 3 ... CpLpCpCiRiCiC2RiC2C3 ■ ■ ■ CpRpC p ciL[ii 

such that the sequences qLiCiL^q, (cjLjCj + i) 2 <j<p and (cjRjCj + i)i<j< p (with c p+ i = C\) code 
disjoint trees, see Figure [2] (left) for an illustration. 

For this type of equivalence class, E[W C ] = r(a 2 ) 1 and there are 

E E (2/i + l)nCat(/,)Cat(r f ) 

p —3 /-| 4-h Ip i=l 

+r 1 +---r v =l-p 

such equivalent classes. Hence, the contribution to the term of order 1 / n in the asymptotic 
expansion of the moment ra 2 / of this subcase is: 

(L E Oh +1) j)cat(/,-)Cat(r,). 

p —3 l\-\ h Ip i —1 

+ri+—r p =l-p 

Finally, let us consider the second case, where edges of the cycle are visited both ways. 
This corresponds to sequences of the form: 

z'iLiCic 2 L 2 c 2 c 3 ... L b _iCb_iCi,L h cijCij-iRi,_iCi,-i ... c 2 cikici 

CpLpCpCp—i... Lh^-iCb+iChRfrChCh+iRb+i... CpRpCpC\L^i\ 
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Figure 2: A graph with one cycle. Left: edges of the cycle are visited one way 
only. Right: edges of the cycle are visited both ways. 


such that the sequences z'lhiCiZ/p'i, (cjLjCj + i) 2 <j< p and (c/fqC/+i)i</<p (with c p+ \ = C\) code 
disjoint trees, see Figure [2] (right). 

For this type of equivalence class, E[W C ] = (o' 2 ) 1 and there are 

E E ptai+nnCatfOCat^) 

p —3 Zi —f-h Ip 1=1 

+ri+---r p =l-p 

such equivalent classes. Hence, the contribution to the term of order 1 / n in the asymptotic 
expansion of the moment m^i of this subcase is: 

l E E p(2'i + l)nCat(/i)Cat(r f ). 

p=3 /jH 1-/p i=l 

+ri+-"r p =Z-p 

Therefore, the total contribution to the term of order 1 / n in the asymptotic expansion of 
the moment mu of this case is: 

m 2\ = \ E (p + r ) E (2/i + 1) n Cat (^) Ca t( r z)- 

p=3 Zq-f (-Zp Z—1 

+ri + ---rp=Z-p 


2.4 Generating series 

We are now interested in the computation of the generating series of the terms of order 1/n 
where the generating series of the Catalan numbers will play a central role. We denote this 
series by 

T( x ) = Cat (k)x k . 

k> 0 
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It has a radius of convergence 1/4 and satisfies the following equalities: 


T(x) = 1 + x(T(x)) 2 = j— 


1 — \/l — 4x 
2 x 


Our computations will also involve the first two derivatives of T, and some basic computa¬ 
tions yield the following useful equations: 


T' = 

jf/ _ 


r 3 


1-xT 2 ' 

IT 5 


+ 


2T 5 


(1 — xT 2 ) 2 (1-xT 2 ) 3 ' 


( 1 ) 

We start with the generating series of mid: 


SiO) = Yj m 2k xk = -y E k ( k + l)Cat(fc)x fc 

fc >0 Z fc >0 

= —^x 2 T"(x) — xT'(x) 
xT 3 

~ “(1 -xT 2 ) 3 ' 


Similarly for m^: 

S 2 (x) = X] = 4 E E ( 2 Pi + 1 ) Cat (Pi) • • • Cat (P4)x k 

k> 0 k>2pi~\ hp4=fc—2 

= V f E (2pi + l)Cat(p 1 )xP i> | f X] Cat(p 2 )x^ 

\P1>0 / \P2>0 / 

= -^x 2 (2xT / + r)r 3 

a x 2 T 5 
a 4 1 — xT 2 ' 

And for 

%(*) = E m 2 k xk = “2 E E ( 2 Pi + 1 )Cat(p 1 )Cat(p 2 )^ 

k >0 k>l pi+p2—k —1 

= ^2 X f E (2pi + l)Cat(pi)xP^ f X] Cat(p 2 )x^ 

u \Pi>Q J \P2>0 J 

s 2 

= -^x(2 xT' + T)T 

_ s 2 xT 3 
cr 2 1 — xT 2 ' 
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And finally for m '^: 


S 4 (x) = Yj m 2 k xk = E E (p + r ) E ( 2J i + 1) ri Cat (^)Cat(r / )x /c 

fc>0 k>3 p=3 /jH-h/p i=l 

+?"iH— r p =k—p 

= E (P + r)xfT 2 r-\2xT + T) = £> + r)xn 2 r 

p>3 L X1 p> 3 

x 3 T 7 x 3 T 7 

~ (l^xT 2 )3 + (2 + r) (l-xT 2 ) 2 ' 

The total generating function of terms of order 1 /ft is thus 


S(x) = S 4 (x) + S 2 (x) + S 3 (x) + S 4 (x) 
xT3 2 „2 t 4 


(x 2 T 4 — 1) + (2 + r) 


x 3 T 7 


+ (2 + r) 


(1 - xT 2 ) 3 
xT 4 

“(1 - xT 2 ) 2 
xT 4 

“(1 - xT 2 ) 2 ' "(1-xT 2 ) 2 
x 3 T 7 xT 3 


(1 — xT 2 ) 2 

x 3 T 7 xT 3 

+ 


+ 


xT 3 


(1 — xT 2 ) 2 1 — xT 2 yu 4 


1 — xT 2 yu 4 
a , S 2 

—rXT + 

cr 1 


(cl 2 S 2 

-^r*E 2 + 


+ 2 - 


x 


3p7 


+ 2 


x 2 T 5 


+ 


xT 3 


+ r 


+ 


(1-xT 2 ) 2 1 — xT 2 y Vu 4 


1 — xT 2 1 — xT 2 
t -2 ) xT 2 + - 1 






We gathered terms in the last equality so that the last four terms vanish in the GUE setting 
for which r = 0, x/V 4 = 2 and s 2 /j 2 = 1. The first four terms therefore correspond to the 
generating series of the terms of order 1 / n for the GUE which we know to be null. It is easy 
to verify that indeed 


xT 4 0 x 3 T 7 n x 2 T 5 xT 3 
(1 - xT 2 ) 2 + (1 - xT 2 ) 2 + 1 - xT 2 + 1 - xT 2 ~~ '' 


and therefore 


S(x) 


= r 


= r 


= r 


X' 


3 j7 


(1 — xT 2 ) 2 
x 3 T 7 


+ 


X 


T 3 


1-xT 2 
xT 3 


«- 2 ) xT 2 + E- 1 


cr 1 


+ 


r ( X T 2 + 1) ^ + 


xT 3 


1 — xT 2 




(1-xT 2 ) 2 1-xT 2 

{T WE_(^ + ( ,r 2 + i)(i-,r 2 )) + r ^ 


-2-r) xT 2 + - 1 — r 




cl „ \ 

—r-2-r) xT 2 + -=■ - 1 - r 


a* 


a^ 


(9) 


Let 

wz)= i T (?H( z ~' /z ^) 

be the Stieltjes transform of the semi-circle law. It satisfies: 

1 -U 2 = U\Jz 2 -4. 
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A simple computation using Q yields 



n n 2 

~ - 4 + Vz 2 - 4 

^ 7 Z' _1_ 

2Vz 2 -4 ^23 

_ r /I / 1 1 

“ 2 U Vz 3 2 + z + 

which is the Stieltjes transform of the measure v in the Theorem. 
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—r -2-r)n +-*-1 -r 


(7 Z 


+ 


w 


CL 


4 J Vz 2 -A V ^ 4 

i \ n 2 
+ 


-r-2-r)H^ + ^-l-r 




\fz 2 — 4 / \/z 2 — 4 


—r — 2 — r) T~L 2 H—— 1 — r 


u A 
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